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SFECTRUMS AND DIFFUSION IN A ROUND TURBULENT JET 
By Stanley Corrsin and Mahinder S. Uberoi 


SUMMARY 


In a round turbulent jet at room temperature, measurement of the 
shear correlation coefficient as a function of frequency (through band- 
pass filters) has given a rather direct verification of Kolmogorof f 1 s 
local-isotropy hypothesis. 

One -dimensional power spectrums of velocity and temperature fluctu- 
ations, measured in unheated and heated jets, respectively, have been 
contrasted. Under the same conditions, the two corresponding transverse 
correlation functions have been measured and compared. 

Finally, measurements have been made of the mean thermal wakes 
behind local (line) heat sources in the unheated turbulent jet, and the 
order of magnitude of the temperature fluctuations has been determined. 


INTRODUCTION 


At the present time there apparently exists no statistical theory 
of turbulent shear flow. One significant theoretical consideration has 
been proposed: The hypothesis of local isotropy, originated by 

Kolmogorof f (references 1 and 2) . Kolmogorof f has suggested that the 
fine structure in turbulent shear flow may be isotropic; recent experi- 
ments of Townsend (references 3 and *0 in a. turbulent wake seem to 
verify this hypothesis. 

In view of this situation, the various experimental researches in 
the field follow two courses: First, they attempt to verify or disprove 

local isotropy; second, they try in all conceivable ways to make measure- 
ments that may shed light upon the basic nature of the turbulent shear 
flow, so that the foundations for a successful theory can be laid. 

During recent years it has become evident that measurements of the 
intensity of turbulence alone cannot provide sufficient information 
about the statistical and dynamical properties of the flow fields. Such 
quantities as correlations, spectrums, probability densities, the various 
terms in the turbulent kinetic-energy balance, and so forth may be 
expected to reveal many essential features of the problem. 
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The types of turbulent shear flow that have come under close 
experimental scrutiny are the boundary layer (references 5 and 6), the 
plane channel (reference 7), the plane wake (references 3 and 4), the 
plane single free-mixing region (reference 8), and the round jet (refer- 
ences 9 and 10) . The present work is a continuation of that reported 
•in references 9 and 10. 

The general objectives of this investigation have been to learn 
something more about the flow in a fully developed round turbulent jet 
and about the heat transfer in such a flow. The work reported here has 
fallen into three phases: (a) An attempt to establish the presence or 

absence of local isotropy, (b) a comparison of velocity- and temperature- 
fluctuation fields when the over-all boundary conditions on mean velocity 
and temperature are effectively the same, and (c) a study of the diffu- 
sion of heat from a local (line) source in the turbulent flow. 

The only specific experimental verification of local isotropy in 
a turbulent shear flow to date was by Townsend in the plane wake behind 
a circular rod ( references 3 and 4) . He found that the skewness and 
flattening factors of the probability density of 8u/8t in the shear 
flow are very closely equal to those in the (effectively isotropic) 
turbulence far behind a grid. Since differentiation emphasizes the 
higher frequencies, his measurement shows, in essence, that the values 
of certain statistical quantities related to the smaller eddies in a 
shear flow are the same as the values for the smaller eddies in isotropic 
turbulence. He also found the microscale of u in the stream direction 
to be nearly \[2 times the microscale of v in that direction, a relation 
which is exactly true for isotropic turbulence. 

Until recently, only mean-velocity and mean-temperature distribu- 
tions were measured to provide a comparison of the transfer rates of 
momentum and of heat in turbulent shear flows with over-all heat 
transfer. The previous report in this round- jet investigation (refer- 
ence 10) included a beginning on the problem of direct comparison of 
the velocity and temperature fluctuations as well as some measurements 
of velocity- temperature correlations. The fluctuations in a warm turbu- 
lent wake have been studied by Townsend. 

Up to the present time, however, there appears still to be no 
successful hypothesis to account for the well-known fact that heat 
(and other scalar quantities, like material) is diffused more rapidly 
than momentum in a turbulent flow. Thus, more detailed study of the 
fluctuations seemed in order. 


■^-Private communications, 1949* 
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The first real analysis on the diffusive property (for scalar 
quantities) of a homogeneous turbulent field was Taylor T s well-known 
work, "Diffusion by Continuous Movements" (reference 11). 

The mean thermal wake behind a line source of heat in a flowing 
isotropic turbulence has been carefully measured by Schubauer (refer- 
ence 12) in the region close to the source, and by Simmons (measure- 
ments reported in reference 13) over an extended range. Taylor *s theory 
of diffusion by continuous movements is directly applicable to the 
diffusion from a line source of heat in a homogeneous turbulent field, 
and he has made a generalization to permit application of the method in 
a decaying isotropic turbulence (reference 13 ). 

In a turbulent shear flow, the only published measurements 6f the 
thermal wake of a local source seem to be those of Skramstad and Schubauer 
in a turbulent boundary layer. 2 These were reported by the experimenters 
(reference lb) and by Dryden (reference 15). The temperature distribution 
across the wake in shear flow is decidedly skew. On the other hand, in 
the isotropic turbulence, it is to all intents and purposes a Gaussian 
curve. Of course, in a homogeneous field this curve (measured close 
enough to the source so that the Lagrangian correlation coefficient is 
still effectively unity) is simply the probability density of the lateral 
velocity fluctuations. This relation may be roughly true for shear flow 
as well, in which case the contrasting results mentioned above would mean 
that the probability density of v(t) is skew in shear flow, but Gaussian 
in decaying isotropic turbulence. 

No measurements had been made of the fluctuations in these thermal 
wakes, and it was felt that some information on the nature of these 
might help further a general understanding of the diffusive process. 

This investigation was conducted at the Aeronautics Department of 
the Johns Hopkins University under the sponsorship and with the finan- 
cial assistance of the National Advisory Committee for Aeronautics. 

The authors would like to acknowledge many stimulating conversations 
with Dr. L. S. G. Kovasznay, and to thank Dr. F. H. Clauser for his 
helpful criticism. Donation of the hot- jet unit by Dr. C. B. Millikan, 
Director of the Guggenheim Aeronautical Laboratory at the California 
Institute of Technology, is greatly appreciated. Mr. Philip Lebowitz 
helped to set up much of the laboratory equipment. 


P 

Wieghardt has recently measured the diffusion of heat from a local 
source in a turbulent boundary layer (reference l6) , but his source was 
flush with the solid surface and thus he was studying a different 
problem, that is, one more directly related to micrometeorological 
conditions . 
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SYMBOLS 

d 

diameter of orifice (l in.) 

X 

axial distance from orifice 

r 

radial distance from jet axis 

u 

axial component of mean velocity 

V 

radial component of mean velocity 

w 

tangential component of mean velocity 

^max 

maximum U at a cross section on axis 

u 

axial component of instantaneous velocity fluctuation 

V 

radial component of instantaneous velocity fluctuation 

w 

tangential component of instantaneous velocity 
fluctuation 

u' S ^ 


V 1 = ft 

• 

e 

instantaneous temperature difference (measured above room 
temperature as reference) 

e 

mean temperature difference (measured above room 
temperature as reference) 

®max 

maximum mean temperature difference at a cross section 


maximum mean temperature in jet at orifice 

•9 

t3 ' = \J~d^ 

instantaneous temperature fluctuation (=6 - 0 ) 

t 

time 

e 

voltage fluctuation 

a, 3 

sensitivity of a diagonal hot-wire to u and to v, 
respectively 
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R uv 

%» v n 

■ u n' = fh? 
v n ' = 

n^uv 

n F u,n F v 


n l 

k 

k l 

k o 

F l( k l) 

F l*( n l) 


shear correlation coefficient (uv/u T v T ) 

instantaneous contributions of u(t) and v(t), 

respectively, in a narrow frequency band of nominal 
frequency n; in a Fourier series discussion, the 
nth harmonic of periodic u(t) and v(t), respectively 


shear correlation coefficient for a narrow band of 

frequencies (u n v n /u n f v n 1 ) 5 this function is referred 
to as the "shear-correlation spectrum" or, briefly, 
the "shear spectrum" 

phase angles 

in Fourier series analysis, n F u = 

/ 00 

n-^v = b n b n ^, where and 

series coefficients 

cyclic frequency in general; in particular, magnitude 
of radial coordinate in three-dimensional frequency 
space 

cyclic frequency of one-dimensional spectrums of u(t) 
and 3(t) 

wave-number magnitude for three-dimensional 
spectrums (2rtn/u) 

wave-number magnitude for one -dimensional 
spectrums ^ 2jtnp fi) 

a reference constant with dimensions of wave numbers 

one -dimensional power spectrum of u(t) in terms of 
wave number 

one -dimensional power spectrum of u(t) in terms of 
frequency 
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F(k) 

Gi(ki) 


Gi»(ni) 


G(k) 


R 


y 


**X 


Ax 


A, 


X 

l 


I 

i 


<t> 

Ti = t/t 
A 


T 


three-dimensional power spectrum of velocity fluctuation 

♦ 

one-dimensional power spectrum of t 9 ( t ) in terms of 
wave number 

one -dimensional power spectrum of $(t) in terms of 
frequency 

three-dimensional power spectrum of temperature 
fluctuation 

transverse correlation function of u measured symmet- 
rically about jet axis 

longitudinal correlation function of u 

transverse correlation function of 3 measured symmet- 




rically about jet axis 

longitudinal scale of u-fluctuations 

longitudinal scale of t 3 -fluctuations 

lateral scale of u-fluctuations 

lateral scale of $ -fluctuations 

longitudinal microscale of u-fluctuations 

lateral microscale of u-fluctuations 

lateral microscale of t 9- fluctuations 

longitudinal microscale of $ -fluctuations 

distance downstream from local heat source in x-direction 

lateral distance, perpendicular to source line, from 
local heat source 

dimensionless temperature ratio ^0/0 max ) 


standard deviation of mean-temperature distribution in 
wake behind local heat source 

pulse spacing 
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h pulse height 

j pulse width 


EQUIPMENT 

Aerodynamic Equipment 


The 1-inch hot- jet unit is shown schematically in figure 1. The 
centrifugal blower is driven by a l/2-horsepower direct- current motor. 
Heat is added through two double banks of coils of No. 1 6 Nichrome wire. 
As can be seen in the sketch, a good part of the heated air is directed 
around the outside of the jet-air pipe in order to maintain a flat 
initial temperature distribution in the jet. A vacuum- cleaner blower 
is used to help the air through this secondary heating annulus, and this 
warm air is fed back into the intake of the main blower. 

The section of relatively high velocity between heaters and final 
pressure box permits adequate mixing behind the grid, to insure thermally 
homogeneous initial jet air. 

Figure 2 is a photograph of the unit as set up previously (refer- 
ence 10); the present arrangement is essentially the same. 

All turbulence measurements were made with an initial jet total 
head in the range from 3.5 to 5.0 inches of water. In free turbulent 
flows there is no detectable effect of jet Reynolds number over a much 
wider range of Reynolds numbers than this. 

When the jet was run unheated, there was a slight temperature rise 
through the blower and duct. In the measurements of thermal wake behind 
a local heat source, correction for this ambient- temperature field was 
unnecessary. For all hot runs, the orifice air temperature was very 
close to 200° C, about 175° above room temperature. 

Three different "local heat sources" were used: 

(a) A straight diametrically strung wire of 0.008-inch Nichrome 

(b) A 2- inch-diameter Nichrome ring 

(c) A 4- inch-diameter Nichrome ring 

Because of the extremely high turbulence levels encountered in a free 
jet, a measurable thermal wake could only be obtained by using source 
temperatures in the range from 300° to 700° C. This undoubtedly led to 
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some local buoyancy effects, but, even with this order of temperatures, 
the thermal wake was barely detectable 1 inch downstream. 

The Reynolds numbers of these heat-source wires were about as 
follows : 

For straight wire on axis, 

150 , based on air temperature 
22 , based on wire temperature 

For 2 -inch-diameter ring, 

110 , based on air temperature 
l 6 , based on wire temperature 

For 4- inch-diameter ring, 

59 y based on air temperature 
9 , based on wire temperature 

No noticeable additional turbulence was generated by these wires, and 
no average momentum defect could be detected with a flattened total- 
head tube, even as close as l/4 inch downstream. 


Measuring Equipment 

The measuring instruments used were: Total -head tube, Chromel- 

Alumel thermocouple, and hot-wire anemometer (also used as resistance 
thermometer) . 

The hot-wires were nominally O.OOO 635 -centime ter platinum, about 
1.5 millimeters in length, etched from Wollaston wire. The etched 
platinum was soft-soldered to the tips of small steel needle supports. 

A discussion of heat loss from a wire at various ambient temperatures 
is given in appendix A. 

The basic hot-wire-anemometry equipment was purchased from 
Mr. Carl L. Thiele of Altadena, California. One of the two identical 
heating circuits is shown in figure 3 * 

The amplifier, with resistance-capacitance compensation network, 
is given in figure 4. The uncompensated gain is constant to within 
±2 percent over a frequency range from 3 to 12,000 cycles per second 
(fig. 5). With the wires and operating conditions used (time constants 
on the order of 1 millisecond) , the over-all compensated response was 
good over the same range. Correct setting of the compensation network 
was determined by superimposing a square wave upon the hot-wire bridge 
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(reference 17). Unfortunately, in a free turbulent shear flow the 
ambient disturbance is so great (because of the extremely high turbu- 
ence levels) that calibration cannot be made in the flow to be studied. 

The vacuum- thermo couple signal output was measured either with a 
millivoltmeter or by the average deflection rate of a fluxmeter. 

The various spectrums reported here were measured with a modified 
General Radio Type 760-A Sound Analyzer (reference 18) . The changes 
in output stage (fig. 6 ) were made to eliminate the direct-current 
component and to obtain linear instead of logarithmic response. As 
modified, the sound analyzer had rather undesirable frequency- response 
characteristics, particularly a day-to-day shift in relative amplifi- 
cation of the higher-frequency ranges. The frequency- response calibra- 
ion in figure 7 is plotted in terms of voltage squared, since this was 
the quantity ultimately measured. 

The frequency pass band for this analyzer is far from the optimum 
rectangular shape. However, the slopes of the two sides are suffi- 
ciently steep that no appreciable error is attributable to noninfinite 
slopes, with the spectrums measured in this investigation. Figure 8 is 
an experimentally determined band shape. There was fair similarity of 
band shape over the entire frequency range. For computational purposes, 
an equivalent rectangular pass band was defined as indicated in the 
figure. 

The instrument is a type recording constant-percent hand width, 
measuring the product of power spectrum times frequency. This has 

obvious advantages in the high-frequency range where there is so little 
turbulent energy. 

Possibly the chief disadvantage of the actual band shape is the 
extremely sharp peak, causing a great deal of fluctuation in the output 
signal, and making any simple meter-reading technique virtually impos- 
sible in the low- and medium-frequency ranges. Consequently an inte- 
grating technique was devised, making use of the negligible restoring- 
torque characteristics of a Sensitive Research Company fluxmeter. The 
integrating technique used is shown schematically in figure 9 . Actually 
a bank of vacuum thermocouples was used, and the resistances shown are 
just typical values. The signal put out by the thermocouples is a 
highly fluctuating direct current. 

The bucking circuit was necessitated by the following combination 
of requirements : 

(a) For the lowest frequencies reasonable consistency could be 
obtained only by integrating over periods as long as 3 minutes 
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(b) Appreciable s tatic bearing friction in the fluxmeter demanded 
more or less continuous motion of the needle 

(c) The restoring torque of the fluxmeter is no longer negligible 
in the range of very large deflection. Hence it was desirable to keep 
total deflection to a minimum. 

Thus, most of the average direct- current component of the thermo- 
couple signal was bucked out, and the constant bucking current was read 
on a precision microammeter- The fluxmeter needle fluctuated more or 
less about the zero-deflection point during the time of integration, 
and its reading at the end of this time ordinarily gave a small correc- 
tion on the result. 

In the highest-frequency range, overload considerations on the 
sound analyzer limited the signal drastically, and only a part of the 
thermocouple direct current was bucked. 

For the determination of average wake temperatures behind the local 
heat sources, the thermocouple voltage was measured with a Leeds & 
Northrup type K-2 potentiometer. 

Oscillograms were taken from a blue oscilloscope tube by means of 
a General Radio Type 65 I-AE camera, using fast film. 


PROCEDURES 
Velocity Spectrum 


The power ( or energy) spectrum of the longitudinal velocity fluctu- 
ations at a point in the unheated jet was measured by conventional hot- 
wire -anemome try technique, with a continuously adjustable band-pass 
filter, as described under EQUIPMENT . 

Temperature Spectrum 

The power spectrum of the temperature fluctuations in the heated 
jet was measured by using the hot-wire effectively as a simple resis- 
tance thermometer (reference 19) . The amplified voltage signal was 
analyzed exactly as in the measurement of velocity spectrums. 
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Shear-Correlation Spectrum 

For the shear-correlation spectrum,, the quantity to he measured is 
the correlation coefficient between a narrow frequency band of 
u-fluctuations and the same narrow frequency band of v-fluctuations, at 
the same point in the flow field. 

The method, for any particular nominal frequency, was to pass the 
various voltage signals (e-^, e 2 , e^ + eg, e^ - eg) from an x-type 

shear- (or v T -) meter through the band-pass filter after amplification. 
By appropriate combination of the mean-square values of these four 
signals (identical, with total-shear measurement), there results 


n R uv = u n v n/ u n'V 


where the subscript n indicates the narrow band of nominal frequency n 
cycles per second. A justification for the validity of this procedure 
is obtainable by considering the two velocity- fluctuation components as 
periodic functions. Of course, this is not a real proof. 

If a symmetrical x-meter is assumed, the two instantaneous voltage 
signals are 


e-^ = au + Pv 

► 

e^ = au - Pv 


(i) 


For periodic fluctuations, 


u = a^ cos (2itnt + 0 n ) 
n=l 


v 


00 


n=l 


► 

b n cos (2rtnt + t n ) 


( 2 ) 


In this simple case, the correlation coefficient for any spectral line 
is merely 


n R U v ~ 


cos (0 n - i|r n ) 


( 3 ) 
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If equations (2) are substituted into equations (l), it is easily- 
shown that 


cos ( 0 n - i n ) = 


n e l " n e 2 


\/n( e l + e 2) 2 n( e l 



(k) 


in complete analogy to the conventional method for measuring R uv with 
an x-meter. The algebraic details are given in appendix B. 


When the meter is not perfectly symmetrical (a^ ^ 09 ; Pi ^ 32 
the formal processes, both algebraic and experimental, become excessively 
involved. Consequently, in actual practice the effect of unavoidable 
unsymmetry in the x-meter was essentially nullified by taking double sets 
of readings at each frequency, rotating the instrument l 80 ° about the 
axis of flow direction between sets. 


Velocity Correlation Function 


The double correlation Ry between longitudinal velocity fluctu- 
ation at pairs of points on opposite sides of the jet axis was measured 
only in the unheated jet. Hence the standard hot-wire -anemometry tech- 
nique was used. The two hot-wires were always equidistant from the 
axis (on a diameter), so that they were under identical operating 
conditions . 


Temperature Correlation Function 

In the hot jet, the wires traversing symmetrically as for Ry were 
operated as simple resistance thermometers, so that the double temperature 
correlation S could be determined directly. 

Mean Temperatures behind Local Source 

To determine mean temperatures behind a local source, traverses 
were made with a Chromel-Alumel thermocouple, whose voltage was measured 
with a Leeds & Northrup type K-2 potentiometer. 

Temperature Fluctuations behind Local Sources 

To determine temperature fluctuations behind local sources, the 
fine platinum wire was operated at small currents, so that it worked 
essentially as a resistance thermometer . 


NACA TN 2124 


13 


EXPERIMENTAL RESULTS 


Mean-velocity and mean-temperature distributions for various 
orifice temperatures are presented in reference 10. 


Shear-Correlation Spectrum 

The spectrum of shear correlation coefficients n R uv = %v^/u n 'v n ' 
was measured in the unheated jet at x/d =20 at a radial station corre- 
sponding to maximum shear at this cross section. Figure 10 shows quite 
definitely that n R uv (n) is a function decreasing monotonically to zero. 
Thus, the hypothesis of local isotropy is seen to be verified in a very 
direct way. The value of the directly measured total-shear correlation 
coefficient R uv = uv/u'v' is indicated in the figure. 


Velocity and Temperature Spectrums 

The one-dimensional power spectrum F 1 (k]_) of the longitudinal 

velocity fluctuations u(t) was measured at two radial positions in 
the unheated jet, at x/d = 20. Figure 11 gives the two spectrums, one 
measured on the axis, and one measured at about the maximum-shear 
location. Plotted against wave number ^k^ = 2itni/uj, the two spectrums 

are identical within the experimental scatter. The solid line drawn as 
approximation to the points is made up as follows: 

(a) For 0 < k^ < 1.25> it is Von Karman's semiempirical formula 
(reference 20): 



Constant 



(5) 


(b) For k^ > 1.25, a nonanalytical curve has been faired in 

The Von Karman expression was used primarily to simplify the 
problem of extrapolation to k-^ = 0 and to shorten the work of trans- 
formation to three-dimensional spectrums. (See "Velocity and Temperature 
Spectrums" under ANALYSIS OF RESULTS.) 

The corresponding one -dimensional power spectrums of temperature 
fluctuations Gp(kx) were measured in the heated jet (0 O = 170° C) 
at x/d = 20. One spectrum is on the jet axis; one is at about the 
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radial station of maximum heat transfer (and shear). From figure 12, it 
can be seen that they differ noticeably in the high-frequency range, but 
are essentially identical in the low- and moderate-frequency ranges. 

The curves used to approximate the experimental points are as follows: 

On the axis - the Von Karman formula is used over the entire range 

At the maximum-heat-transfer point - 

(a) For 0 < <1.0, the Von Karmen formula is used 

(b) For kj >1.0, a nonanalytical curve has been faired in 

Figures 13 and 14 contrast velocity spectrums with temperature 
spectrums at corresponding radial stations. 


Transverse Correlation Functions 

The double correlation function Ry = ^l^ 2 / u i ' u 2 T y measured symmet- 
rically about the axis at x/d = 20 in the unheated jet, is plotted 
in figure 15- Of course, since the wires are in identical flow condi- 
tions, u^ T = U 2 * = u ! (say), and Ry = u^^/u^. 

The double correlation function Sy = * 8 1*8 2 /® 1*“® 2* ^he heated 

jet at x/d = 20 is plott ed in figure 1 6 . Since this was also measured 
symmetrically, Sy = * 8 3 ^ 2 /^* where 8p T = * 82 1 = 1 (say). 

Clearly, the range of measurable temperature correlation exceeds 
the range of measurable velocity correlation by an amount greater than 
can be attributed simply to the fact that the hot jet is wider than the 
unheated jet (reference 10) . 


Mean Thermal Wakes behind Local Heat Sources 

Typical radial distributions of average temperature behind a 
straight (diametrical) wire, a 2- inch-diameter ring, and a 4- inch- 
diameter ring at x/d = 20 in the unheated jet are shown in figures 17, 
l8, and 19, respectively. The points in these figures are not direct 
experimental points, but merely serve to indicat^ the faired results 
for different distances downstream. All these have been corrected for 
the previously mentioned small ambient-temperature field in the unheated 
jet. There was rather large scatter (illustrated only in fig. 20) due 
to the small temperature differences measured and to the extremely large 
degree of fluctuation present. The 4- inch ring is slightly outside of 
the fully turbulent jet core (reference 9); consequently the results for 
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this case are not of direct interest in a study of fully developed 
turbulence. The figures show that each of the thermal wakes possesses 
similarity well within the accuracy of measurement. 

All of the thermal wakes spread linearly in the measured range 
(figs. 21, 22, and 23) • From Taylor's theory of diffusion by continuous 
movements, this simply indicates that, for the maximum downstream station 
studied, the Lagrangian correlation coefficient of the v-fluctuations 
has still not departed appreciably from unity. 

For a straight-line source at the jet axis, because of conservation 
of heat, it follows immediately from similarity and linear spread that 
the maximum temperature at a cross section in the wake 0 max must 

decrease hyperbolically with increasing downstream distance. 

Let £/^ max = $(?])> where t \ = £/& and 6 is some characteristic 

width of the wake, for example, the £ at which 0 = l/2; conservation 
of heat gives 

00 

6 d£ = Constant (6) 


where the mean- velocity changes are neglected. Then, 

^max 5 = Consta nt/li (7) 

where 

noo 

J l = dT l 

Jo 

The same is true of the annular wake in the range of I so small 
that 5 « r. 

Single traverses were also made behind a straight-line heat source 
for two other cases: 

(a) With the line source set perpendicular to r at a radius of 
1 inch, a temperature traverse was made in the r-direction, 

at | = l/2 inch (fig. 24) 

(b) With the line source set on a diametral line, a temperature 
traverse was made perpendicular to r at a radius of 1 inch; 

£ = l/2 inch (fig. 25) 
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Temperature Fluctuations behind Local Heat Source 

As can be anticipated, the temperature fluctuations close behind 
a local heat source are quite different in nature from the velocity 
fluctuations at the same point or from the temperature and velocity 
fluctuations in a turbulent flow with over-all heat transfer. Since a 
suitable source produces no additional turbulence, 3 the velocity fluc- 
tuations should be the same as in undisturbed flow. On the other hand, 
the " turbulent” thermal wake close to the source must be simply a very 
narrow laminar thermal wake which is fluctuating in direction as v(t) 
fluctuates. Since all of the fluid outside of this unsteady laminar 
wake is of constant temperature and the temperature fluctuations can 
only be positive, the general character of the oscillogram of ( t) in 
figure 2 6 is understandable. These records were taken about 3/8 inch 
dojmstream from the straight-line heat source, and about 3/l6 inch off 
the wake axis. All of the oscillograms were made with insufficient 
compensation for the hot-wire thermal lag, in order to suppress the 
/high frequency) noise and thus permit the basic form of T9(t) to 
stand out. From these two oscillograms of u(t) and check measurements 
of the turbulence levels for the two cases, it appears that the source 
wire has made no appreciable change in the turbulence. 

Measurements of the intensity of the temperature fluctuations 
across a section at | = 0.4 inch are given in figure 27- In the same 
vicinity the values of u’/U and v'/u are on the order of 20 percent. 
The extremely high values of i3 T /F are n °t surprising since the mean 
temperature difference is due only to the presence of the fluctuation. 
Since the simple -re si stance -thermometer theory is based upon the assump- 
tion of small fluctuations compared with absolute temperature (refer- 
ence 19 ) y it is expected that these measurements are about as accurate 
as the measurements of much lower t9 ! /0 in the hot jet (reference 10) . 


ANALYSIS OF RESULTS 
Shear-Correlation Spectrum 


A very convenient check upon the shear-spectrum measurements can 
be gotten by the simple expedient of computing the total (or ’’net”) 
turbulent shear correlation coefficient (which was also directly 
measured) from this spectrum and the turbulent-energy spectrum. Again, 
an elementary Fourier series treatment serves to justify (not prove) 

3An ideal source would also produce no average momentum defect. 
However, as mentioned previously, the momentum wakes of the local sources 
were relatively so small as to be completely undetectable as close as 
l/k inch downstream. 
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the intuitive idea that the total correlation coefficient R, 


'UV 


IS 


simply a weighted average of the 


n^uv 


(i.e., the cosines of the phase 


angles), weighted simply hy the product of the square roots of the 
energy spectrums of u and v. The calculation is given in appendix B, 
and yields the relation 


R 


■uv 


00 - / 

= ^ (n F u n^v) 1 2 n^uv 


(8) 


n=l 


Unfortunately, 25 cycles per second is the lower limit of the 
measured frequency range, so that some extrapolation must be made to 
lower frequencies which contain much of the turbulent energy. Since no 
theoretical basis yet exists to guide this extrapolation (like the 
Von Karman formula in the case of energy spectrums), guesses had to be 
made as to the maximum and minimum of reasonable-looking extrapolation 
curves. These are plotted in figure 28, along with the energy spectrum 
of the u-fluctuations . Since no spectrum of the v-fluctuations was 
measured, the expression actually used for computing is 

*uv = / F i( k i) k^uv a*i (9) 

Jo 

The best of the three extrapolations tried (extrapolation Q)) 
gives R uv = 0.46, which is satisfactorily close to the directly 
measured value of 0.44, especially since* there is no reason to suppose 
that the spectrum of the v-fluctuations is identical with the spectrum 
of the u-fluctuations, except in the range of local isotropy. 


Velocity and Temperature Spectrums 

The one-dimensional spectrums of velocity and temperature fluctu- 
ations, as plotted in figures 11 and 12, respectively, are area- 
normalized; that is, they are defined such that 


F l( k l) = 1*° 


( 10 ) 


G 1 (k 1 ) dkx = 1.0 


( 11 ) 
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However, the original measurements were made on an absolute -value 
basis, so that the total fluctuation levels u'/u and d'/© could 
be used as checks on the spectrums . The spectrums as measured 
were Fp*(ni) and Gi*(np), defined such that, ideally, 


Fi*( n i) dnp = u 2 


( 12 ) 


G l*( n l) 


JO 



( 13 ) 


Integration of F-j* as indicated in equation (12) yielded the 
following: 


On the axis (r = 0) , 




1/2 


F^* dn^ 


= 0.28 


(u'/U = 0.22, directly measured) 


In the maximum-shear region (r = 1.0 cm), 

1/2 


II 


F^* dnp 


= 0.52 


(u'/u = O.lO, directly measured) 


Similar integration of the measured temperature spectrum in the 
heated jet yielded the following: 


On the axis (r = 0) , 


1/2 


fI 


G ± * dn 1 


= 0.21 


(d'/© = 0.18, directly measured) 
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In the maximum-shear region (r = 4.8 cm). 


1 

6 




0 . 4l 


(t3 '/© = O. 36 , directly measured) 


On the jet axis, -where conventional small -perturbation hot-wire 
theory may still be moderately accurate, the agreement is satisfactory. 

It_should be noted in passing that these directly measured values 
of - 3'/ 6 are appreciably higher than those reported in reference 10 . 

No explanation for this difference is apparent. 

The longitudinal scale of u-fluctuations is obtainable approxi- 
mately from the power spectrum: 



(14) 


which follows from the fundamental Fourier transformation, 


F i( k i) ■ t 


R x cos ( kpx ) dx 


Jo 


(15) 


in the limit k^_ — > 0 . 

An analogous treatment of a temperature -fluctuation field leads to 
an identical expression for the longitudinal scale of d-fluctuations : 


A x = f G X ( 0 ) ( 16 ) 

Within the accuracy of measurement, the longitudinal scale of u(t) 
was found to he the same on the jet axis and in the maximum-shear 
region: 


L x = 3*6 centimeters 
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The same was true of the- longitudinal temperature scales in the hot jet: 

A x = 2.4 centimeters 


In all these computations the measured spectrums were extrapolated 
to zero wave number with a parabola, as first suggested by Dryden 
(reference 21) . Approximately the same numerical values are obtained 
with the Von Karman formula illustrated in the plotted curves. 


In order to compare L x in the unheated jet with A x in the 
heated jet, L x may be multiplied by 1.15, the jet-width ratio 
at x/d = 20 for these two initial temperatures (reference 10) . The 
n corrected M L x is then 4.15 centimeters. 


It must be recalled that the Fourier transformation relation 
between time spectrum and space correlation would be exactly true only 
if the turbulent fluctuations at a point were due to pure rectilinear 
translation (by U) of a fixed fluctuation pattern. For the free jet 
flow, the extremely high turbulence levels make such a transformation 
very uncertain. Therefore, equations (l4) and (l6) (and the resulting 
scales) can only be considered as crude approximations. 

The longitudinal microscale, 



( 17 ) 


may also be computed approximately from the one -dimensional power 
spectrum, by 




ki 2 Fi(ki) dk x 


(18) 


(Primes signify differentiation when applied to correlations and spectrums.) 
Again, an analogous approach to temperature fluctuations gives 



1 

2 


U 


0 


k l 2c l( k l) <^1 


( 19 ) 
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Unfortunately, the high-frequency ranges of the Gp ' s are too 

uncertain to permit reasonable extrapolation and the use of equation (l 9 ), 
although.it can be seen from figure 13 that A. x > Z x . However, equa- 
tion (l 8 ) has been used to compute the longitudinal velocity microscale. 
Figure 29 is a plot of the integrand of equation (l 8 ) . The integration 
and appropriate computation give 

\ x = 0.44 centimeter 

and the same value for both radial positions. 

Since the turbulence on the axis of such a jet seems to be rather 
isotropic (the experimental evidence is that uv = 0 and u' ~ v'), 
the lateral microscale 


\ = 


2 

Ry"(0) 



( 20 ) 


is of the order of that is, 

\ £ 0.31 centimeter 


With the assumption of isotropic turbulence on the jet axis, it is 
possible to compute the three-dimensional power spectrum F(k) from 
the one-dimensional spectrum F^ki). Heisenberg (reference 22 ) has 
given the inverse transformation 



and the desired F(Fp) is readily found to be 

F(ki) = 2k 1 [kiF 1 "(k 1 ) - F^ki)] (22) 


A three-dimensional spectrum computed in this way is given in figure 30. 
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The corresponding spectral transformations for the three-dimensional 
isotropic fluctuation field of a scalar quantity (temperature, for 
example) are simply (reference 23 ) 

^00 

G l( kl ) = § nr (23) 

J k l 

and 

G(k x ) = (24) 

With the assumption of isotropic temperature fluctuations on the 
jet axis, equation (24) has been used to compute G(k) . 

The three-dimensional velocity and temperature spectrums on the 
jet axis are very nearly the same - provided that the isotropy assump- 
tion is reasonably good. The curve in figure 30, which shows the 
general nature of both F and G, is simply the transformation of the 
Von Karman approximations to F^_ and G^_. 

It must be recalled, however, that F^ is in the unheated jet, 

while Gq_ is in the heated jet. Although the results of references 9 
and 10 indicate no essential change in the detailed dynamics of jet 
turbulence as a result of moderate increase in jet temperature, there 
is an appreciable increase in jet width at a given x/d. As mentioned 
earlier, the width ratio between 6 Q = 175° and 6 Q = 0° is about 1.15 
at x/d = 20. 


Transverse Correlation Functions 

The integral area under the velocity correlation function (fig. 15) 
may be considered to give a sort of lateral scale of turbulence in the 
jet, although the result is not associated with any particular region 
in the jet. The conventional expression, 




00 


gives a scale, Ly 


0.67 centimeter 


Ry(r) dr 


(25) 
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If a lateral temperature scale is defined in similar fashion, 


A y = 


S y (r) dr 


(26) 


then for the heated jet at x/d = 20 it turns out that Ay = 0.77 centi- 
meter from the function as given in figure l6. 


Appropriate comparison of these transverse scales may he had if Ly 
is multiplied by the jet-width ratio: l.l5Ly = 0.77 centimeter, the 

same value as Ay. However, the two correlation functions that yield 

these net areas are still quite different in shape. The contrast is 
shown in figure 31* Clearly, even though the net areas are identical, 
there is nonzero temperature correlation over appreciably greater 
distances . 


A rough approximation to the microscale of turbulence can be gotten 
by guessing at the osculating parabola for Ar = 0. In this particular 
case "guessing” is more appropriate than "fitting," since the job is 
entirely extrapolatory in nature. Figure 32(a) shows the vertex region 
of Ry with the parabola that corresponds to 


X = 



= 0.28 centimeter 


(27) 


This value is in surprisingly good agreement with the 0.31 centimeter 
obtained from the power spectrum on the axis. In fact, the agreement 
must be regarded as fortuitous, since the difference is appreciably less 
than the experimental uncertainty. 

If the temperature-fluctuation field is again considered analogously, 
the transverse microscale of temperature fluctuations (fig. 32(b)) is 


Z 


2 

Sy"(0) 


1/2 


= 0.4-3 centimeter 


(28) 


For comparison, 1.15*. = 0.32. 


Mean Thermal Wakes behind Local Heat Sources 

The rate of spread of the thermal wake close behind a local source 
of heat was first used by Schubauer (reference 12) as a means of 
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measuring the intensity of lateral velocity fluctuations v'/u. A 
detailed discussion of this technique has been given by Taylor (refer- 
ence 13 ) and need not be repeated here. The results of such a computa- 
tion, compared with direct x-meter measurements of v'/u, are as 
follows : 


(in.) 

(=) 

' U / wake 

(=) 

/x-meter 

0 

0.18 

0.185 

1 

.22 

.20 

2 

.40 

.30 


Th e x-meter measurements were corrected for the effects of both u' 
and uv upon the slightly unsymmetrical meter. 

It is possible to get some additional information about the fluctua- 
tion field by computation from the turbulent-heat-transfer equation. In 
particular, an estimate of the distribution $v across a section of the 
thermal wake behind the line source may be made as follows. 


The steady turbulent-heat-transfer equation for low velocity (negli- 
gible viscous dissipation to heat) , negligible molecular heat conduction, 
and constant density is, in Cartesian tensor notation. 


u, 


1 chc^ 



U 9 ) 


For the region in the immediate vicinity of the jet axis, assume 
that condi tion£ approximate those in a homogeneous field of turbulence 
that is, V = W = 0 and U = Constant = U max> 

Then equation ( 29 ) becomes simply 

1 * ■ fi r si) - ! 


(30) 
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The assumption of small turbulence level implies -3u« 9U max and 

hence 


u ^ 

u max 



(31) 


This would be a good approximation in turbulence far behind a grid 
placed in a uniform stream, hut is certainly rather crude here. 

The final assumption, that of similarity in the thermal wake, is 
well supported by the experimental results. Then let 

r~ - *(& s <32) 

°max ^ 

where A is the standard deviation of the mean-temperature distribution, 
and, according to the theory of diffusion by continuous movements, is 
therefore proportional to the standard deviation of the probability 
density of v(t) as well. Specifically, for a small | in a homo- 
geneous turbulent flow, 



(33) 


Similarity also implies that 


t9v - ^xnax^max (say) (3^) 


Equation (7) may be written 


0 max A = Constant 


(7a) 


Then, with equations (32), ( 3^-) > and (7a), equation (3l) maybe trans- 
formed to 


*5 =. ^ A. (,f) 
dT| a | dt| 


or, with equation ( 33 ) > 


doj 

dt] 


v' 

'"max 



(35) 


( 36 ) 
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Since V/U max is constant in this approximation, 


0>(t)) 


U, 


max 


f(t)) + Constant 


and the boundary condition, cjo = 0 at T) = 0, gives finally 

fly v 1 £ 6 

^max^max U max A 6 mav 

which is conveniently written in the form 


(37) 


dv £ 

e u el 

max max max 


(38) 


In figure 33 this function is plotted against £/| for the traverse 
l/2 inch downstream from the straight-wire heat source, across the jet 
center. 


DISCUSSION 
Local Isotropy 


The monotonic decrease to zero in shear correlation coefficient 
with increasing frequency seems to he decisive evidence for the exis- 
tence of local isotropy at sufficiently high Reynolds numbers. It is 
interesting to note that the spectral region of negligible shear 
(n^ > 1000 cps in the present particular determination, for example) 
contains only about 1.5 percent of the turbulent kinetic energy in u2. 
Of course, this is by no means an indication of the importance of the 
existence of local isotropy in a turbulent shear flow. A more perti- 
nent comparison would be with figure 29, which shows in effect dissi- 
pation as a function of frequency. From this it appears that about 
90 percent of the dissipation of turbulent kinetic energy to heat takes 
place in essentially isotropic turbulence. This permits the use of the 
Taylor expression for dissipation in isotropic turbulence (reference 24). 
Of course it also implies that the isotropic relation between longi- 
tudinal and lateral microscales, \ x = ^A,, will be fairly accurate 

even in the region of high turbulent shear. Furthermore, it implies a 
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universal dimensionless spectral function for all turbulent flows, in 
the high-frequency region. For turbulence at this Reynolds number, it 
appears that a universal part of the spectrum exists only for k]_ > 7-9> 
which is well beyond the point of slope -^. In fact, the ^ point" 

in this spectrum is just at k-^ = 1.0. 


At this point a few remarks on the appropriate type of measurement 
for verification of local isotropy may be in order. In particular, a 
careful distinction must be made between the shear spectrum as 

presented in this report and the power spectrum of the randomly fluc- 
tuating quantity uv which might be measured with a multiplying circuit 
followed by a frequency analyzer. 

Local isotropy specifies that restriction to a sufficiently small 
domain in a turbulent shear flow shows up isotropy in the various statis- 
tical properties that are studied within that domain. It implies that 
all three components of the frequency (or wave-number) vector must be 
large, and not merely the magnitude of the vector. Clearly then a good 
indication of isotropy is zero correlation between orthogonal velocity- 
fluctuation components; this means that the highest-frequency parts 
of u, v, and w are uncorrelated with each other. Hence it is clear 

that if u n Vn decreases to zero, with increasing frequency, faster than 
the product u n f v n ' decreases to zero, local isotropy exists. In terms 
of coefficient, this merely requires that n ^uv decrease to zero 
eventually. 


Now consider the fluctuating quantity uv. In a turbulent shear 
flow uv f 0, so that uv consists of a direct-current component with 
superimposed random fluctuations. Since the conventional electronic 
techniques eliminate the direct current, the quantity to be analyzed 
would be uv - uv as a function of time. If local isotropy were 
present, the lower frequencies of u and v would be rectified in the 
multiplying process, and therefore the oscillogram and power spectrum 
of uv - uv would have relatively great emphasis on the high frequencies. 
In other words, if the naively measured power spectrum of uv were used 
as an indication of local isotropy, it would show a trend opposite to 
that of u n v n ; that is, it would decrease more slowly than the 

product In general, the measurement of n Ru V seems like a 

much more specific and direct approach than the measurement of the 
power spectrum of uv. Presumably, a (somewhat more complicated) 

Fourier series discussion like that in appendix B could also be carried 
out for the power spectrum of uv. 
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Velocity and Temperature Spectrums 

The apparent identity of' the velocity power spectrums Fp(k]_) on 
the jet axis and in the region of maximum shear is only approximate and 
has been determined only down to kp = 0.1. There still exists the 
possibility of measurable divergence in the lowest wave-number range. 

The good degree of agreement indicates that, in diffusing from the region 
of maximum production (near the maximum-shear region) to the region of 
maximum dissipation (on the jet axis), the turbulent kinetic energy has 
not done any gross migrating in the wave -number space. 

On the other hand, the apparent decided difference between temper- 
ature power spectrums measured on the axis and in the maximum-heat- 
transfer region seems to indicate such a migration. However, the 
considerable scatter at the highest measured frequencies renders 
definite conclusions impossible. 

Somewhat more specific conclusions can be drawn from the comparison 
between one -dimensional velocity and temperature spectrums. On the jet 
axis, for example, in spite of distinct differences between these two 
spectrums, it turns out that within the experimental scatter (which is 
considerable) the three-dimensional power spectrums may be much more 
nearly identical. The fact that they did in fact come out to be identi- 
cal over a wide range of wave number when computed' from the empirically 
fitted Von Karman formula must certainly be regarded as pure chance. 

This is true not only because of the experimental uncertainty, but also 
because these spectrums were measured in two similar but different flows, 
whose characteristic lengths probably differed by 15 percent. 


Kinematic and Thermal Scales 

From the extrapolated zero -wave -number intercepts of the one- 
dimensional spectrums, the following longitudinal scales were obtained 
at x/d = 20: 

L x = 4.15 centimeters 
A x = 2.4 centimeters 


This L x is 15 percent greater than the unheated-jet value, to allow 
for the greater width of the heated jet. Thus, L x /A x =1.7- In a 

homogeneous, isotropic field of velocity and temperature fluctuations, 
it turns out (reference 23) that, if the three-dimensional power 
spectrums of velocity and temperature are proportional, L x /A x = I.50. 
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It may also be noted that, if the measured ratio were in an isotropic 
field, Ly = |L X and A y = A x , so that L y/ A y = O. 85 . The ideal value 

would be 0.75- Actually, the integrals of the transverse correlation 
functions R y and S y are considerably less than the scales that 
would be expected, according to these relations, in a homogeneous 
isotropic turbulence. 

On the other hand, the relative values of longitudinal and lateral 
kinematic microscales follow the isotropic relation, \ x = \[2\, at 
least within the experimental uncertainty. This is on the order of 
±25 percent in the case of the parabola "fitted" at the vertex of R y . 

Unfortunately, the temperature spectrum on the jet axis is not 
extended sufficiently far to permit computation of longitudinal micro- 
scale Z x there. The spectrums in figure 13 show only that Z x is 

considerably less than l x ; that is, Z x is considerably less than 
0.44 centimeter. It may be remarked in passing that isotropy for a 
scalar quantity means equality of longitudinal and lateral correlation 
functions. The lateral microscale, Z = 0.43 centimeter, obtained by 
fitting a parabola at the vertex of S y seems of reasonable magni- 
tude relative to 1.15X. =0.32 centimeter. 


Transverse Correlation Functions 

Of course, the reason L y and A y as determined by integration 
of functions R y and S y are not related isotropically to L x and A x 

is that over most of the range of Ar the probes are in decidedly non- 
isotropic turbulence. Thus, there is no reason to expect L v = 1 l v 

or A y = A x , when L x and A x are computed from the spectrums. 

An examination of the behavior of these two symmetrically measured 
correlation functions shows that there is nonzero correlation over a 
considerable part of the jet, but that the relatively small scales result 
from the rather extensive regions of negative correlation. This behavior 
is emphasized by a comparison of R y with the corresponding function in 
some typical isotropic turbulence downstream of 1-inch-mesh grid (refer- 
ence 25) • Figure 34 shows the contrast clearly. 

It is conceivable that such an extended region of negative corre- 
lation is characteristic of turbulent shear flow. However, until some- 
one establishes this in a shear flow whose transverse extent is very 
large compared with the maximum correlation distance, it may be safer 
to guess that the "excess" amount of negative correlation is simply due 
to a slight irregular waving of the jet as a whole. In reference 9 it 
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was assumed that, since Ry actually goes to zero at large values 
of Ar, there is no over-all "whipping" of the jet. However, such a 
conclusion does not appear to he completely warranted. 

In the section entitled "Transverse Correlation Functions" under 
ANALYSIS OF RESULTS, it was found that in the heated jet Ly ~ Ay On 

the other hand, it is well known that the lateral rate of transfer of 
heat is appreciably greater than the lateral rate of transfer of 
momentum, as was first found by Ruden (reference 26) from mean-velocity 
and mean -temperature measurements. Since diffusion is essentially 
Lagrangian in nature, while L and A are Eulerian scales, the above 
results are not necessarily in contradiction. The appreciably greater 
distance over which Sy ^ 0 (as contrasted with Ry) may, however, be 

related to the fact that the mean thermal jet diameter is appreciably 
greater than the mean momentum jet diameter. 


The mean-temperature distribution close behind the straight-line 
heat source on the jet axis is effectively symmetrical, and closely 
resembles a Gaussian curve in shape (fig. 20) ; this shows that the 
probability density of v(t) on the axis is more or less Gaussian, as 
in isotropic turbulence. 

The me an -temperature distributions close behind the two ring heat 
sources are decidedly skew. However, some of this skewness seems to 
be due simply to the curvature of the line sources. Therefore, the 
temperature distribution across the wake of a straight wire set tangent 
to the circle r = 1 inch was measured. Neglecting the effects of 
mean-velocity gradient, this curve (fig. 24) is proportional to the 
probability density of the radial velocity fluctuation v(t) in the 
shear region. It is seen to be slightly skew; the skewness factor 


is computed directly from this curve. The thermal wake measurements of 
Skramstad and Schubauer behind a line source in a turbulent boundary 
layer (reported in reference 15) show a skewness of 0.38* The dif- 
ferences in sign and magnitude of these two skewness factors suggest 
lateral turbulence -level gradient as the cause. The gradients in v'l U 
are of opposite sign in these two flows. 


Probability Density of v(t) and w(t) 
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Calculation from figure 25 shows that the probability density of the 
tangential fluctuation w(t) is symmetrical. It may be noted that on the 
axis of such an axially symmetric flow there is no distinction between 
radial and tangential velocity fluctuation; hence figure 20 also applies 
to w(t) on the axis. 


Temperature Fluctuations behind Local Heat Source 

The extremely high temperature -fluctuation levels ( fl'/ET > 1.0) 
encountered in the wake of the line heat source are easily understood 
from a brief consideration of the nature of the temperature field. 

Close behind the source, there is just a single narrow laminar thermal 
wake which is being blown in random deviations from the ^-direction by 
the turbulent fluctuations. The gross turbulent thermal wake is simply 
the wedge-shaped region over which this relatively narrow wake wanders. 
Hence the total thermal signal at any fixed point in the gross wake 
consists simply of a series of pulses, where each pulse corresponds to 
an occasion upon which the laminar wake swept over the point. Obviously, 
the frequency of occurrence of pulses will decrease monotonically with 
increasing transverse distance from the center of the gross wake. 

If this type of temperature signal is represented schematically by 
periodic square pulses of height h, width j, and fundamental wave 
length t (fig. 35) > then it can be easily deduced that the fluctuation 
level is 



where it is recalled that t9 = G - 6 by definition. 

Two pulse spacings of interest are 

(1) t = 2j; then -d'/G = 1.0 

(2) t— > 00 ; then -d '/o’ — > °° 

Hence, the measured results for $ '/ ~6 seem quite reasonable in 
both order of magnitude and in qualitative behavior across the gross 
thermal wake of the local heat source . 

The distribution of ^v/^max^max computed from the thermal wake, 
plus the availability of the measurements of d T /^max an( ^ v V^max' 

suggests the computation of the heat-transfer correlation coeffi- 
cient tSv/tS'v*. Unfortunately, when the results of figures 33 and 27 
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are used, a part of the correlation-coefficient distribution reaches 
impossible values (slightly above unity).. It must be concluded that 
the absolute values of $v / GmaxUmax are too inaccurate for such a 
computation. 

Sources of Error 

Aside from the specific instances mentioned earlier in this 
section, the sources of experimental error are much the same as outlined 
on pages 27 to 28 of reference 10. Additional uncertainties arise in 
the spectrum measurements, especially in the higher-frequency range, 
because of (a) rapid changes in the calibration of the sound analyzer 
(band peak response against frequency, fig. 7) and (b) slight static 
friction of fluxmeter bearings. 

In general, it should be emphasized that measurements by conven- 
tional ( small -perturbation) hot-wire anemometry in a flow of this high 
level of turbulence cannot be considered as accurate absolute -value 
measurements. Even on the jet axis, where the level is a minimum and 
conditions are relatively steady, there is no reason to believe that 
absolute values are better than within, say, ±10 percent of the "correct” 
values. However, relative behaviors are undoubtedly determined, and 
dimensionless measures of the type of correlation coefficients are more 
accurate than absolute values. 

None of the measurements reported here have been corrected for 
finite length of hot-wires. 

SUMMARY OF RESULTS 


From measurements in a round turbulent jet at room temperature of 
the shear correlation coefficient as a function of frequency, of 
velocity and temperature fluctuations with and without jet heating, 
and of the mean thermal wakes behind local heat sources, the following 
statements may be made: 

1. The Kolmogoroff hypothesis of local isotropy is verified for 
the shear flow in a round turbulent jet. This is concluded from the 
monotonic decrease to zero of the shear-correlation spectrum (u n v n /u n ' v n ’) 
with increasing frequency n. 

2. The one-dimensional power spectrums of longitudinal velocity 
fluctuations and of temperature fluctuations appear to be basically 
different. 
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3- The three-dimensional power spectrums of velocity and temper- 
ature fluctuations on the jet axis seem to he roughly alike - if the 
assumption of isotropy in this region he true. It may then follow that 
the difference in the one -dimensional power spectrums is a direct 
manifestation of the fact that velocity and heat are vector and scalar 
quantities , respectively. 

4. The ratio of longitudinal to lateral scale (for hoth velocity 
and temperature fluctuations) is considerably larger than would follow 
from isotropy. Longitudinal scales are measured on the jet axis, while 
lateral scales involve a traverse of most of the fully turbulent core 
of the jet. 

% 

5. The ratio of longitudinal to lateral kinematic microscale on the 
jet axis is about equal to the isotropic value. 

6. The longitudinal thermal microscale (from one -dimensional power 
spectrum s) is less than the longitudinal kinematic microscale, but the 
lateral microscales (from correlation measurements) have the opposite 
relation; that is, the thermal is greater than the kinematic. 

7. The probability density of the radial fluctuation v(t) on the 
jet axis is effectively Gaussian. The probability density in the shear 
region is slightly skew. 

8. The temperature-fluctuation field in the wake behind a local 
heat source consists of a randomly waving narrow laminar thermal wake. 
Hence the temperature signal at a fixed point is a random-pulse type 
of function. Its fluctuation intensity is on the order of 100 percent 
on the center line, and increases toward the edges. 


The Johns Hopkins University 

Baltimore, Md., August 17, 1949 
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APPENDIX A 


HEAT LOSS FROM A WIRE AT VARIOUS AMBIENT TEMPERATURES 


In figure 3 of reference 19 a rough check vas made on the 
temperature-variation first term in King's (reference 27) equation 
for the steady static heat loss from a cylinder perpendicular to a 
fluid stream, at low Reynolds numbers. The conventional form is 


where 


and 

R 

R a 

% 

l* 

d* 

cc* 

k* 

C P 

P 

i 


i 2 R r= 

= A + BwU 

R - R a 


(Al) 


A 


Z*k* 

Roa* 


B 


c 2 


l* 

RqCC* 



wire resistance 

wire resistance at ambient fluid temperature 
wire resistance at 0° C 
wire length 
wire diameter 

temperature coefficient of change of resistivity of wire 
material 

thermal conductivity of fluid at ambient temperature 
specific heat of fluid at ambient temperature 
density of fluid at ambient temperature 
current 


c l» c 2 


empirical constants 
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In reference 19 the check on A as a function of temperature was 
made by assuming the second terra in equation (Al) to be exact in its 
temperature variation. Then each measured calibration point at any 
velocity and temperature led to a value for A. 

The present check was carried out more completely; a full calibra- 
tion curve was run for each ambient temperature. From this, both A 
and B were determined. Figure 3 6 gives the results compared with 
King's predicted variation, using physical constants from reference 28. 
Each point corresponds to a calibration. The vertical line through a 
point obviously does not represent the over-all uncertainty; it simply 
shows the range of values that could be gotten by drawing different 
reasonable-looking straight lines through the same set of original 
calibration points. From the figure it can be seen that King's equation 
predicts the temperature variation of A quite well. The changes in B 

( the slope of the calibration line in the plot of -i— — against 
\ ^ “ ^a J 

are so small that the experimental scatter is as great as the changes 
predicted for these temperature differences. 
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APPENDIX B 

MEASUKEMENT OF SHEAR -CORRELATION SPECTRUM 


The two voltage signals from an ideal symmetrical x-meter are 


e-^ = ou + 3v 
eg = otu - 3v 


(Bl) 


Suppose that the velocity fluctuations are periodic: 


u - X a n cos (2jrnt + 0 n ) 
n=l 

00 

v = ^ h n cos (2jmt + \|/ n ) 
n=l 


(B2) 


Of course, there would he no loss in generality if 0 n or \|/ n were 
taken as zero. 


The quantity to be measured is 


R UnVn 

= u n' v n' 


(B3) 


For two simple harmonic functions the correlation coefficient is sinrply 
the phase angle. Thus, 


n R U v - cos (0 n - ^ n ) 

Substitution of equations (B2) into equations (Bl), followed by 
trigonometric transformation, gives 

oo_ 

e i = 

¥ cos 0 n + 3b n cos ) cos 2irnt - 

(aajj sin 0 n + pb n sin \|r n ) sin 2rtntj 


(B4) 


(B5a) 
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00 


e 2 - ^ ^(^n cos 0n - P b n cos ^n) cos 2nnt 

- 

( <xa n 0n ~ 3^n ^n) sin 2irntJ 

(B5b) 

When these two signals are put separately through a narrow band- 
pass filter that passes only the nth harmonic, the two output voltages 
may be represented as 

n e l = *[(0*41 cos 0n + P b n cos *n) cos 2irnt - 


(oca^ sin 0 n + 3b n sin i n ) sin 2*nt] 

(B6a) 

n e 2 = K^o^n cos 0n - P b n cos \4 n ) cos 2imt - 


(oa^ sin 0 n - |3b n sin t n ) sin 2«nt] 

(B6b) 

where K is an attenuation factor. 


For brevity, write 


n e l = K ( A n cos 2*nt - B n sin 2jrnt) 

(B6c) 

n e 2 = K(c n cos 2nnt - D n sin 2itnt) 

(B6d) 


These filtered, signals go next into the vacuum— thermocouple unit, 
which puts out the mean-square values. 


~ n e l 2 = -^n 2 cos y (2irnt) - 2A n B n cos (2irnt) sin (2rtnt) + 


B n ^ sin^ (2jrnt) 
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1- n e 2 ^ = C n 2 cos 2 (2irnt) - 2C n D n cos (2jrnt) sin (27rnt) + 


D n 2 sin 2 (2imt) 


where K ! is an over-all attenuation factor. 


But, cos^ x sin^ £ and cos sin ^ 0, over a large number of 


wavelengths . Thus , 


n e l 2 K ^ ( A n 2 + Bn 2 ) 
Ifi? * ^ ( c n 2 * D n 2 ) 


(B7) 


Then, within the approximation, 


e 2 e 2 _ a 2 p 2 , ^ 2 D 2^ 

n e l n e 2 ” 2 \ n Un + u n ) 


(b8) 


and when the expressions for A, B, C, and D are substituted, it 
turns out that 


n e l 2 " n e 2 2 = 2K 1 a0a n b n cos (0 n - ^ n ) 


(B9) 


The necessity of determining a and 3 is ordinarily avoided with 
a symmetrical meter, if only the correlation coefficient is required. 

The sum and difference of the two wire voltages are 


e l + e 2 = 2au = 2a / a n cos (2jtnt + 0 n ) 

1 


e l - e 2 = 23v = 23 > h n cos (2jmt + \|r n ) 

1 


(BIO) 
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Filtering gives 


n( e l + e 2) = 2K aa n cos (2nnt + 0 n ) 
n( e l “ e 2) = 2KP4) n cos (2imt + i n ) 
Passage through the vacuum thermocouple gives 

n( e l + e 2) 2 = 2K'a 2 a n 2 


(Bll) 


n( e l - e 2 ) 2 = 2K'P 2 b n 2 

Combination of equations (Bll), (B9), and (B4) gives the final 
result: 


n^uv - 


~ 2 p 2 
n e l " n e 2 


n( e l + e 2) 2 n( e l - e 2) 2 ] l/2 


(B12) 


The computation of total-shear correlation coefficient from shear- 
coefficient spectrum suggests itself as a useful check possibility: 


R uv = uv/u'v* 


(B13) 


vith the Fourier series for u and v, 


u' 2 = u2 = iXI *n 2 


(B14) 


.,2 


oo 
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and the instantaneous cross product can he transformed to 


00 oo 


UV = II a^bju^cos 0 n cos 2imt - sin 0 n sin 2«nt) x 

n=l m=l 


(cos cos 2«nt - sin t m sin 2«mt) 


The time average of this expression is 

anM 0013 cos ^n + sin 0n sin ^n ) 


— 1 
UV = — 

2 n=l 


or 




Vn cos ( ^n “ *n ) 


(B15) 


Thus, 


u'v' R uy - 2_ Bn* 3 n^uv 

1 


In terms of the Fourier coefficients, 

00 

y ^t^n n^u 

*\iv = 


— 2 V~ 2 \ 1/2 


(B16) 


But n F u = a n y / a^ and n F y = b n 2 j / By, 2 are simply the normalized 
one-dimensional energy spectrums of u and v, respectively. Therefore, 


hiv 


n=l 


/ \ l/2 
(n F u n F v) n^uv 


(B17) 
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Figure 2.- The jet unit. 
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Square -wave generator input circuit 
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Figure 3.- Control circuits 
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Figure 4.- Hot -wire -anemometer amplifier. 
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Figure 5.- Amplifier frequency response without compensation. 
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(b) Modified output circuit. 

Modification of General Radio Type 760 -A Sound Analyzer. 


Figure 6.- 


Output 

(Output at 60 cps) 
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Figure 7.- Response of sound analyzer with modified output circuit. 


NACA TN 2124 


v_n 

ro 



Figure 8.- Filtration characteristics of sound analyser with modified output circuit. Area of rectangle 

equals area under curve. 
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Figure 10.- Variation of shear correlation coefficient n R uv with frequency. Round turbulent jet. x/d = 20; 

U = 270 centimeters per second. 
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Figure 11.- One -dimensional power spectrums of u(t) measured in 1-inch 
unheated jet at x/d = 20. Computed scales: L = 3.3 centimeters for 

-X 

both stations, A x = 0.31 centimeter for both stations. 
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Figure 12.- One -dimensional power spectrums of •a(t) measured in 1-inch 

heated jet at x/d = 20. Computed scales: A x = |g 1 (0) = 2.18 centimeters 
for both stations. 
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Figure 13.- One -dimensional spectrums of u(t) and i3(t) on axis of 1-inch 

jet at x/d = 20. 
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Figure 14.- One -dimensional power spectrums of u(t) and s(t) in maximum - 

shear region of 1-inch jet at x/d = 20. 
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Figure 15.- Symmetric transverse correlation of 
x/d = 20 in 1-inch unheated jet. 


u, measured about axis at 
R y = 



Figure 16.- Symmetric transverse correlation of measured at x/d = 20 in 1 -inch heated jet. 

s y = vi/^. 
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Figure 17.- Similarity of temperature behind line source of heat on jet axis at x/d = 20 in unheated jet. 

?! = ? at which I - |5 mSK . 



Figure 18.- Similarity of temperature distribution behind 2-inch-diameter-ring source of heat in unheated 

jet at x/d = 20. ^ = ? at which 0 = |e max . 
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Figure 19.- Similarity of temperature behind 4 -inch-diameter-ring source of heat in unheated jet at 

x/d = 20. t lt 5 2 = 5 at which e = |e max - 


ON 

UO 


M VOVM 


Figure 20.- 
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Experimental scatter; temperature behind line source of heat on jet axis, t = 1.27 centimeters; 

x/d = 20. 
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Figure 21.- Spread of heat from a line source in center of unheated jet. x/d = 20. 
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Figure 23.- Spread of heat from 4 -inch-diameter-ring source in unheated jet. x/d =20. = 5 at 


which 9 = - 0 

2 max 
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Figure 24.- Temperature behind line source of heat in unheated jet. x/d = 20; 5 = 1/2 inch. Line source 
set perpendicular to r at a radius of 1 inch. Traverse made in r-direction. 
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Figure 25.- Temperature behind line source of heat in unheated jet. x/d = 20; 5 = 1/2 inch. Line source 
set on a diametral line. Traverse made perpendicular to r at a radius of 1 inch. 
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Figure 26.- Oscillograms of velocity and temperature fluctuations. 
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Figure 27.- Temperature fluctuations behind line source of heat in center of unheated jet. x/d = 20; 

I = 0.4 inch. 
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Figure 28.- Shear spectrum and power spectrum in 1-inch unheated jet at x/d = 20 in region of 

maximum shear. 
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Figure 29.- Viscous dissipation and microscale A on axis of 1-inch jet. x/d = 20; = j k-^F^QcOdk-].; 


A = 0.31 centimeter. 
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Figure 30.- General nature of three-dimensional spectrums on jet axis, x/d = 20. 
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Figure 31.- Comparison of transverse correlation functions. 1 -inch heated jet. 
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(a) Lateral microscale of u-fluctuations. (b) Lateral microscale of -fluctuations. 

Figure 32.- Estimates of microscales. 1-inch jet. x/d = 20. 
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Figure 33.- Thermal wake behind straight-line heat source. Temperature -velocity correlation computed from 

®Amax* 5 = °- 5 inch - 
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Figure 34.- Transverse velocity correlation functions. Ry. = u-^Ug/u^ in round jet; g = ^Ug/u^ in 

isotropic turbulence. 
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Simulation of temperature signal close behind local heat source. 


Figure 35.- 
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Figure 36.- Variation of hot-wire constants with air temperature. 

( ) r , room temperature. 
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